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If neutrinos are Majorana fermions, the lepton mixing parameter space consists of six mixing 
parameters: three mixing angles and three CP-odd phases. A related issue concerns the physical 
range of the mixing parameters. What values should these take so that all physically distinguishable 
mixing scenarios are realized? We present a detailed discussion of the lepton mixing parameter space 
in the case of two and three active neutrinos, and in the case of three active and N sterile neutrinos. 
We emphasize that this question, which has been a source of confusion even among "neutrino" 
physicists, is connected to an unambiguous definition of the neutrino mass eigenstates. We find that 
all Majorana phases can always be constrained to lie between and n, and that all mixing angles 
can be chosen positive and at most less than or equal to tt/2 provided the Dirac phases are allowed 
to vary between — -k and n. We illustrate our results with several examples. Finally, we point out 
that, in the case of new flavor-changing neutrino interactions, the lepton mixing parameter space 
may need to be enlarged. We properly qualify this statement, and offer concrete examples. 



I. INTRODUCTION 

It is now established that, similar to quarks, leptons mix The phenomenon of fermion mixing was identified 
long ago (see, for example, 0, Q). It is a consequence of the fact that we have (at least) three sets of each "type" 
of chiral fermion field (Ql,Lli&r,ur,. . . ) and, while gauge interactions cannot tell different generations apart, the 
interactions responsible for fermion masses can. 

In the quark sector of the standard model, it is well-known that the ambiguity in defining the different quark chiral 
fields allows one to describe the entire flavor sector of the quarks in terms of 10 parameters: six quark masses (three 
for the up- type quarks, three for the down- type quarks) and four real parameters in the Cabibbo-Kobayashi-Maskawa 
matrix (say, the Wolfenstein parameters A = sxnOc, where 9c is the Cabibbo angle, A, p, rj). It is also well-known 
(but often not emphasized) that one can choose all quark masses positive and all quark mixing angles to lie in the 
first quadrant (for the Wolfenstein parameters, this translates into A > 0, A > 0). 

In the lepton sector the situation is similar, with a couple of "twists:" the neutrinos may be Majorana fermions, 
and the neutrino sector of the Standard Model can be easily enlarged to included so-called sterile neutrinos v s . Unlike 
fourth-generation quarks or charged leptons, sterile neutrinos can arise from Standard-Model singlet chiral fermion 
fields, and mix with the so-called active neutrinos after electroweak symmetry is broken. Not only are sterile neutrinos 
of virtually all masses allowed by all current neutrino data, they are natural consequences of several well-motivated 
scenarios that explain the origin of the Majorana neutrino masses, including most versions of the seesaw mechanism 
It has recently been emphasized that there are plausible reasons to expect the "seesaw" sterile neutrinos to be 
low-mass, propagating fields [5|. 

In the absence of sterile neutrinos, it is well-known 0, 0] that the flavor sector of the leptons is described in terms 
of 12 parameters: three charged-lepton masses, three neutrino masses, and six real parameters in the lepton mixing 
matrix, often referred to as the Pontecorvo-Maki-Nakagawa-Sakata matrix. The last are parameterized in terms of 
three real mixing angles (#12, #23, #13) and three complex phases. If the neutrinos are Dirac fermions, two of the three 
complex phases are unphysical, while if one of the neutrino masses is identically zero, one of the three complex phases 
can be "rotated away." 

In this paper, we discuss the physical range of the lepton mixing parameters (mixing angles and complex phases) 
assuming that the neutrinos are Majorana fermions and allowing for the presence of sterile neutrinos. By "physical 
range" wc refer to the following question: 'given a proper definition of the different neutrino states, what values should 
the mixing parameters take in order to describe all physically distinguishable values of all observables?'. Subsets of 
this discussion, usually concentrat ing on the mixing angles and the "Dirac" phase, can be found elsewhere (see, for 
example, [!, d, S EH EH El) EE EH EH for an incomplete list of references) , and many of the results we present 
here in detail can be found throughout the literature. It is, however, our experience that the topic of the "physical" 
lepton-mixing parameter space is often source of confusion and misunderstanding, even among seasoned neutrino 
aficionados. A coarse scan of the literature reveals several incorrect (but, fortunately, mostly harmless) statements 
and misunderstandings. For this reason alone, we think such a discussion is useful. Furthermore, a detailed discussion 
of the physical parameter space for Majorana phases and the physical range for all parameters when there are sterile 
neutrinos is, to the best of our knowledge, new and, we hope, useful. 

Several subtle aspects of the definition of fermion mixing never surfaced with great prominence in the quark sector, 
because (a) both quark flavor and mass eigenstates can be produced or detected, and (b) quark mixing angles turn 
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out to be very small. In the lepton sector, even if one ignores the possibility of Majorana neutrinos or sterile neutrino 
states, the situation is different. For all practical purposes, in the laboratory, all neutrinos are produced as flavor 
eigenstates, and all neutrinos are detected as flavor eigenstates. Moreover, lepton mixing angles are large (or even 
maximal) and preempt the most straight forward definition of the different neutrino mass eigenstates. 

This paper is organized as follows. In Sec. |TT1 we carefully define the mixing matrix and describe how unphysical 
field redefinitions define an equivalence class of physical mixing matrices. We then discuss in detail the two, three and 
many Majorana neutrino cases. We illustrate some of our results with well-known (and not so well-known) examples. 
In Sec. IIII1 we discuss how the situation would change in the presence of new neutrino interactions. We discuss 
both lepton-number violating and lepton-number conserving new neutrino interactions. In Sec. IIV1 we summarize 
our results and offer some concluding remarks. Three appendices contain, respectively, a detailed discussion of the 
notation we use to describe unitary matrices of different dimensionality, an example of a "non-standard" parameter- 
space choice that may be useful in the case of Dirac CP-conservation, and approximate expressions for neutrino mixing 
matrices in the case of three active plus one or two sterile neutrinos. 



II. PARAMETERS AND PHYSICAL RANGES 



After electroweak symmetry breaking, assuming three sequential generations of leptons that couple to the Standard 
Model gauge fields, the Standard Model Lagrangian augmented by Majorana neutrino masses, H v SMi contains 
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Here, £ a (e Q ), a = e,fj,,T, are the left-handed (right-handed), charged fermion fields, i — 1, 2, . . . , N + 3 are 
left-handed mass-eigenstate neutrino fields with a well-defined mass, and v a , a = e, fj,, r, s\, . . . , sjv are left-handed 
flavor-eigenstate neutrino fields that couple diagonally to the left-handed charged leptons. N = 0, 1, . . . parameterizes 
the number of sterile neutrinos. Note that for the charged-current and neutral-current neutrino interactions, the sums 
are restricted only to the so-called active neutrino flavors, v e ,v lM ,v T 0]. 

We choose to work on the weak basis where the charged lepton and neutrino masses (m e ,m fi ,m T and 
mi, m2, ■ ■ ■ |7773 + 7v, respectively) are real and positive. Neutrino mass eigenstates and flavor eigenstates are related 
via the neutrino mixing matrix U: 
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where U is a unitary 3 + N x 3 + N matrix: UdUif — S a /3, U a iU*j = Sy. The indices i, j = 1,2, .. .3 + N and 
a, (3 = e, fx, t, si, . . . , sjv- We further choose to parameterize the lepton mixing matrix as (see [HI, E3] for a more 
detailed discussion) 
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where P'(</>) is a diagonal 3 + N x 3 + N matrix whose diagonal entries are all unity except for the i, i element, given 
by e 4 * (see Appendix [5J. U' is a unitary matrix that results when one chooses a parameterization such that as many 
diagonal complex phases are defined to "to the left" and "to the right" of U as possible.* U' is non-diagonal and 
depends on the henceforth defined mixing angles 9 and Dirac phases S. Diagonal redefinitions of the l a and e a fields 
that preserve m a real and positive render all <fi a values physically equivalent^ The same is not true of the henceforth 
defined Majorana phases <pi (note that we have already "spent" an overall neutrino field redefinition to render the 
(1,1) element of the matrix on the right-hand side of Eq. (|II.3|) equal to unity). The only remaining field redefinition 
that leaves rrii real and positive is to change the sign of v% fields: Vi — ► —v*. Hence, without loss of generality, all U 



will be parameterized as U 



u'(o,g) (n 



3+N ■ 



* If neutrinos were Dirac fermions, U' would contain all of the physically meaningful mixing parameters. 

t There is one subtlety here. In the case of N ^ 0, there aren't any sterile charged-fermion fields around to field-redefine some of the 
</> a . Nonetheless, it is easy to see that different values of these phases do not lead to physically distinguishable observables as these 
parameters do not appear in either the charged-current or neutral-current neutrino interactions and are hence unobservable. 
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We haven't spent all of our ability to redefine neutral and charged lepton fields. From the discussion above, two 
matrices U(6, 6, <pi) and 11(9', 8', ^) describe the same phenomena if 

/3+N \ /3+N \ 

u{e',8',<j>' i )= y[[ p'co.ttJJ \ J[ p j (o,7r)J , (n.4) 

where the diagonal matrices on the left- and right-hand sides have entries that are either 1 or —1. Eq. (|II.4[) will be 
used to constrain the physical parameter space of the mixing angles 9, the Dirac phases S and Majorana phases <pi, 
once the neutrino mass eigenstates are unambiguously defined. 



A. The Two Neutrino Case 



If there were only two active neutrinos (for concreteness, 24 and v^), U would be a 2 x 2 matrix. According to our 
conventions, U is parameterized by 

° = ( Z Z ) = R12 "" p2< * ) - ( Zl ) ( i S* ) • < IL5 » 

where the matrix R 12 is defined in Appendix [A] 

The physical parameter space for 9 and <f> depends on the definition of neutrino mass eigenstates and neutrino 
flavor eigenstates (or weak eigenstates). The definition of neutrino flavor eigenstates is unambiguous: the electron 
neutrino v e is the chiral field that couples to the electron and the PF-boson fields, while the electron is the lightest 
charged- fermion. The same goes for v^. It is the neutrino that couples to the muon, which in turn is the heaviest 
charged-lepton. The neutrino mass eigenstates v\ and need to be similarly defined. Here, we discuss two different 
definitions. Note that both uniquely specify the neutrino mass eigenstates. 

A2. v\ is the lightest state, while is the heaviest state: m\ < 771.2. 

B2. v-i is the state with the largest v e content, while v 2 is the state with the largest i/„ content. The v a content of 



is define to be \U ai \ 2 : \U el \ 2 > \U,~ 12 



Direct examination of Eq. (|II.5|1 reveals that, at most, 9 e [— n, w[ and </> e [—tt,tt[ (case A2) and 9 E [— tt, — 3tt/4] U 
[-71-/4, tt/4] U [3tt/4, tt[ and (j> <E [-7T, tt[ (case B2). 

Eq. pi.4[) is now used to further constrain the physical ranges of 9 and <\>. This is done by identifying, given a fixed 
value of 9 and 0, all 9' and 4>' that satisfy 

R 12 (0')P 2 (0') = N f R 12 (0)P 2 (</))N m , or (II.6) 
R 12 (6»')P 2 (</>') N m = N f R 12 (6»)P 2 (</)). (II.7) 

where Nf = I, P 1 (7r), P 2 (7r), —I and N m = I,P 2 (7r), where I is the identify matrix (2 x 2 in this case). The results 
of all these equations (some of which are redundant) are easy to obtain using the algebra presented in Appendix [Al 
and are tabulated in Table [U For example, 

-IR 12 (6»)P 2 (0) = R 12 (0 + tt)P 2 (</>) = R 12 (6» + 7r)P 2 (</>)I, (II.8) 

so that N f = -I and N m = I lead to 0' = 9 + n and <p' = 4>. 

The table (and the example above) reveals that 9 and 9 + tt are equivalent. Hence, it is enough to restrict oneself 
to 9 € [— 7r/2,7r/2] (case A2) or [— 7r/4,7r/4] (case B2). Furthermore, 9 and —9 are also equivalent. This allows one 
to further restrict 9 £ [0,7r/2] (case A2) or [0,7r/4] (case B2). Finally, the Majorana phase is independently invariant 
under <fi — > <fi + tt. It can therefore be constrained to (f> £ [0, 7r] (cases A2 and B2). This result agrees with a recent 
analysis of leptonic rephasing invariants [l5| . 

The reason for the different physical parameter space for 9 in case A2 and B2 has nothing to do with the redundancies 
of the fields and parameters, but is a consequence of our definition of the mass eigenstates. It is clear that both choices 
describe all different physics scenarios equivalently well. This is accomplished by noting that, in case B2, one must 
allow for both mi < 7712 and mi > mi. This is not true in case A2. It possible to map case A2 into case B2. One can 
write 

U(0, 0) = P 2 (7r)R 12 (7r/2 - 0)P 2 (-0)R 12 (7r/2)e 4 *, (II.9) 



or 




= , e** cos(7r/2-tf) sin ( 7r/2 - 6) ]{ , " — (II . 10 ) 
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TABLE I: Values of 0' and <j>' that satisfy Eq. |jL7t , for different values of N f and N„ 
are the field redefinitions corresponding to the different choices of Nf and N m . 



(see text for details). Also tabulated 
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-I, I 

-i, p 2 W 
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1/2 — » —V2 

f — > — / 

■ ~£p, V2 -> 



-I/ 2 



-lu, V2 



f, + 7T 
-6>, + 7T 

-6, <$> 
+ TY,(f> 
l +7T, -0 



Hence, a mixing matrix with < 9 < tt/A can be mapped into another mixing matrix with 7r/4 < 9 < tt/2 if the 
sign of cf> is reversed and if one "exchanges" v\ <-> i>i (up to an allowed redefinition of the left-handed lepton fields) . 
This means that points (9, ra\,m-z, 4>) in the case A2 parameter space associated to 7r/4 < 9 < tt/2 arc equivalent to 
the points (tt/2 — 9, rri2, mi, —<fi) in the case B2 parameter space. This will be concretely illustrated in the examples 
below. 

The probability P™ c of a neutrino produced as a v e with energy E to be detected as after propagating a distance 
L in vacuum is 

F e 7 = sin 2 20sin 2 (^), (11.11) 

where Am 2 = m\ — m\. In case A2, Am 2 is positive-definite, while in case B2 it can be either positive or negative. 
It is clear that this observable does not depend on the Majorana phase (j>. It is also clear that 9 and —9 and it + 9 
lead to the same oscillation probability. Furthermore, in case A2, 9 and tt/2 — 9 also yield the same value of P™ c for 
fixed Am 2 . This implies that in case B2, for fixed 9, Am 2 and —Am 2 should also yield the same value of P™° . This 
is, of course, the case. 

The probability p™ attcr of a neutrino produced as a v e with energy E to be detected as after propagating a 
distance L in a medium characterized by a constant electron background is 

sin 2 29 . , / Am 2 L 



p ™"" ^ nrv^SB si " 2 ^Vi + *-Mco.»j , ,ii,2, 

where A — 2y/2EGrN e / Am 2 , Gp is the Fermi constant, and N e is the local electron number density. As before, 
Am 2 = to 2 — m 2 . As in the case of Eq. (jll.llj) . the Majorana phase <f> plays no role in determining P™ attcr , and one 
can explicitly see that 9 and —9 and 9 and tt + 9 lead to the same value of p™ attor . i n case A2, 9 and 7r/2 — 9 lead 
to different values of p™ atter for fixed Am 2 since 9 — > tt/2 — 9 leads to cos 29 — > — cos 29 and p™ atter depends on the 
sign of the cos 29 term Q. In case B2, one can see that Am 2 and — Am 2 also lead to different values for p^ attcr fo r 
fixed 9. Furthermore, as advertised earlier, the points (6, Am 2 ) and (tt/2 — 9, Am 2 ) in case A2 are equivalent to the 
points (9, Am 2 ) and (9, —Am 2 ) in case B2. 

The value of the Majorana phase will only affect processes where lepton number is violated. As an example, consider 
the rate for neutrino-antincutrino oscillations, related to the probability that a neutrino produced as a v e with energy 
E is detected, after propagating some distance L, as what appears to be a z7 M . It is [l7l Il8| 

sin 2 29 ( o 2 ~ ( Am 2 L 



Pep oc ^ 2 jmf + m l 2 - 2m x m 2 cos I — — 2<j) j j . (11.13) 

As in the previous cases, it is clear that 9 £ [0,7r/2] covers all distinct values of this neutrino-antineutrino oscillation 
probability since P e p, is only a function of sin 2 29. One can also see that <f> and tt + </> lead to the same expression. In 
the case A2 parameterization, one sees that for fixed m\ and m 2 , 9 and tt/2 — 9 lead to the same expression. This 
is also true using parameterization B2 and exchanging mi mi and <fi — * —<p. To the best of our knowledge, this 
"invariance" was first noticed in 1171. 



5 



B. The Three Neutrino Case 



In the case of three active neutrinos and no sterile neutrinos (N 
conventions, U is parameterized by 



0) U is a 3 x 3 matrix. According to our 
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where, to preserve space, we take advantage of the shorthand notation cy = cos( 



i j and Sjj 



sin%, ij = 12,13,23. 



As in the previous subsection, the question of the physical parameter space for the mixing angles #12, $13, #23 and 
the CP- violating phases <5, 02,03 depends on the proper definition of the neutrino mass eigenstates, v\, ^2, and u 3 . 
Here we discuss three cases: 

A3. v\ is the lightest state, 1/2 is the second lightest state, 1/3 is the heaviest state: m\ < TO2 < m 3 . 

B3. v\ is the state with the largest v e content, V2 is the state with the second largest v e content, and v 3 is the state 



with the smallest v e content: J7 e i > \U e 2\ > \U, 



e:3 



C3. The v\ and 1/2 states are defined such that v\ is always lighter than V2 and ro 2 . — m\ is, in magnitude, the 
smallest mass-squared difference. In this case, mi < 1712 < rn 3 if Ato 2 3 > and m 3 < mi < m,2 if Am^ < 0. 
The former is referred to as the 'normal' mass hierarchy, while the latter is the 'inverted' mass hierarchy. Note 



that we define Am 2 = 



9 2 



Direct examination of Eq. (|II.15j) reveals that, at most, 9ij E [— tt,tt[ (ij = 12,13,23) and 02,03,<5 € [— tt, tt[ 
(cases A3 and C3). For case B3, the situation is more complicated. Naively, #23 € [— 7r,7r[ (ij — 12,13,23) and 
02,03,i$ S [— 7r, 7r[, while sin 2 0i2 < 1/2 and tan 2 6*13 < sin 2 $12. Furthermore, since sin0i3 and 8 always appear as 
the combination sin0i3e ±l<5 we are allowed to either constrain sin #13 > or 8 >G [0,7r]. One can also check that 
Eq. (|II.15[) is unchanged if one simultaneously redefine all mixing angles — * fty + tt, and 5 — » 8 + tt. Both of these 
redundancies are included in Table [TT] (see below). 

In order to further constrain the mixing parameter space, we search for all 9^,8' , 2 , 4>' 3 (ij = 12, 13, 23) that satisfy 



[/(^.,02,0 3 ,<5') = N f t/(%,0 2 ,0 3 ,<5)N m! or 
[/(%0 2 ,0 3 ,<5')N m - N f t/(%,0 2 ,0 3 ,<5)- 



(11.16) 
(11.17) 



Here, N f is equal to I, P^tt) (i = 1,2,3), R«(tt) = P^P^tt) (ij = 12,13,23), or -I = P 1 ( 7 r)P 2 (7r)P 3 ( 7 r). On 
the other hand, N m is equal to I, P*(7r) (i — 2,3), or R 23 (7r) = P 2 (7r)P 3 (tt) . Before proceeding, we need to qualify 
the types of solutions we are looking for. We are only interested in linear solutions of the type = $ + k, where 
k is constant and $ is a generic mixing parameter. These solutions allow one to relate disconnected regions of the 
parameter space and can be used to reduce the volume of the original space. We do not rule out the possibility that 
other solutions exists, but these will only apply to "lower-dimensional" subspaces of the parameter space and are not 
useful when it comes to defining the physically distinguisable values of the mixing parameters. 
For example, 



U(dij,8, 02, <fe)R 23 M = [R 2 

= t* 1 
= [R 2 



h 3 )P 3 (-8/2)P 1 (8/2)R 13 (e 13 )F 3 (8/2)P 1 (-8/2)R 

,12//i \t>2/ J, 1 — \ I>3 



12/ 



' 12 )P 2 (0 2 )P 3 (03)] R 23 W, 

13 )P 3 (<5/2)P 1 (-<5/2)R 12 (^i2)P 2 (02 + ^)P 3 (0 3 + tt)] , (11.18) 

23 + ^)P 3 (-<5/2)P 1 (^/2)R 13 (-0i3)P 3 W2)P 1 (-^/2)R 12 (-0 12 )P 2 (02)P 3 (0 3 )] , 

(11.19) 



where we made use of the notation and the product rules described in Appendix [Al Eqs. (|II.18III.19[) indicate that 
the mixing matrices parameterized by {612, #13, #23, 02, 03}, {#12, #13, #23, <5, 02 + tt, 03 + ?>"}> and {—#12, -#13> #23 + 
tt, 8, 02, 03} describe the same physics. Many other solutions for 9^,8', 2 , 3 are tabulated in Table ITT1 including all 
of the independent ones. Several more, none of which provide any extra constraints, have been omitted. 

Table |TTJ reveals that for fixed mixing angles and 8, matrices with 02 and 02 + tt are physically indistinguishable. 
This also holds for 03 and 03 + tt. Hence both can be constrained: 02, 03 G [0, 7r] (all cases A3, B3, and C3). Keeping 
all parameters fixed, one can easily show that matrices with 612 and 612 + tt are identical, the same holding for 023 
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TABLE II: Values of 

^12,13,23 an d 8, 4>2, 03 that satisfy Eq. (JTTTTtJi , for different values of Nf and N m (see text for details). Also 
tabulated are the field redefinitions corresponding to the different choices of Nf and N m . 
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and 023 + tt and for #13 and — (#13 + tt) . This allows one to further constrain the ranges for the mixing angles to Oij G 
[-71-/2, tt/2[ (ij = 12,13,23) (cases A3 and C3) and 9 12 € [-tt/4, tt/4], 9 13 € [- arctan(| sin^i 2 |), arctan(| sin6> 12 |)], 
^23 € [— 7r/2,7r/2[ (case B3). Finally, points with 9^,5 and —9^,5 + n (for any particular ij = 12,13,23) are also 
physically indistinguishable. We will then choose all mixing angles positive, which implies that the physical range for 
the Dirac phase is 6 € [— tt, 7t[. We are left with 

• cases A3 and C3: 9\ 2 , #13, #23 G [0,7r/2], 5 € [-7r,7r[, 02,03 G [0,7r]; 

• case B3: 6>i 2 G [0,7r/4], 6»i 3 G [0, arctan(sin^i 2 )], 6» 23 G [0,7r/2], (5 g [— 7r, 7r[, 02,03 G [0,7r]. 

This choice is not unique, and is a consequence of allowing S to take all of its naively allowed values. It also agrees 
with the choice made in virtually all of the neutrino literature and seems to be the most natural one. We find one 
circumstance where a different choice may be advantageous from a pragmatic point of view. Using the parameterization 
above, the case of Dirac CP-conservation is contained in two "disjoint" subsets of the parameter space: S = and 
5 = tt. On the other hand, if one chose #13 such that negative and positive values were included, S could be constrained 
to lie within [— 7r/2, 7t/2]. In this case, the Dirac CP-conserving subspace is continuous and, e.g., in case C3, #13 would 
lie within {—it/2, n/2]. This observation was first made by the authors of [HI]. A concrete numerical example of this 
is depicted in Appendix [Bj 

As in the two-neutrino case, one can map one definition of neutrino mass eigenstates into another. Case C3 turns 
out to be the one that is widely utilized in the literature [l|, with case A3 a distant second Unlike the two 

neutrino case, the map between different cases is more involved. Here, for illustrative purposes, we discuss briefly how 
one can relate case A3 and case C3. 

If the neutrino mass hierarchy is normal (in case C3 language, Amf 3 > 0), case A3 and case C3 are identical. If the 
mass hierarchy is inverted (in case C3 language, Am^ 3 < 0), however, the situation is more interesting. Fig. Q] depicts 
two identical copies of a particular point in the lepton parameter space (charged- fermion masses not included), one 
labeled according to the conventions of case A3 (right-hand side) , the other according to the conventions of case C3 
(left-hand side). In order to convert from case C3 — ► case A3, one need only "relabel" the mass eigenstates as follows: 
3 -> 1, 1 -»• 2, 2 -»■ 3. 

Such a relabeling will be accompanied by a non-trivial mapping of mixing angles. Fig. [T] reveals, for example, that 
^e3 12 ^ = ^ii 23 \ which implies that s 1 g 12 ''e l ' 5C ' = 0^3 1 where the parameter superscripts indicate the mass 
labeling (312 is an inverted mass ordering in case C3, 123 is the same picture in case A3). We proceed, using the 



t Case B3 was inspired by a comment made by Alexei Smirnov to one the authors. To the best of our knowledge, case B3 has not been 
explicitly considered in the literature. We can't, however, rule out the possibility that it has been implicitly assumed in several occasions. 
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FIG. 1: Identical copies of a particular point in the lepton parameter space (charged- fermion masses not included), one labeled 
according to the conventions of case A3 (right-hand side), the other according to the conventions of case C3 (left-hand side). 
See text for details. 



language of Appendix \Kl to formally complete the mapping. Specifically, we study the effect of 
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on the mixing matrix. Acting on the mixing matrix, this has the following effect: 
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(11.21) 



The entire transformation may be absorbed into shifts of the physical mixing parameters (plus re-identifications of the 
Majorana phases) up to an overall field redefinition of all charged-lepton fields. Formally, the transformations shown 
in Eqs. (|II.21[) are very simple and can be summarized as a relabeling of the mixing planes. This leads to a situation 
where the complex rotation governed by #13 and S now acts in the 1 — 2 plane, while #i 2 governs the 1 — 3 rotation 
— the "roles" of 12 and l3 are exchanged. To express the result of Eq. (|II.21[) in the form of Eq. 15[) . however, 
one must "commute" the 1 — 2 and 1 — 3 rotations and properly redefine the # 23 mixing angle. We find the following 
mapping of mixing parameters from the 312 (Fig. [ljleft), case C3) to the 123 (Fig. [TJ right), case A3) pictures. 



(312) 

Am 2 3 
Ato| 3 
Ato 2 2 
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-03, 



(11.22) 



where capital-letter parameters (Ci 2 , A, $ 3 , etc) stand for parameters in the 312 picture (Fig. [ljleft), case C3) while 
lower-case letter parameters (ci 2 ,(5, 3 , etc) stand for parameters in the 123 picture (Fig. [IJ right), case A3). These 
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transformations are unique up to trivial sign redefinitions, provided the overall nonphysical rephasing by e l ^ + ^ 3 \ 
The transformations are relatively simple for the 9\ 2 and #13 parameters, but the redefinition is less attractive for 
$23 due to its induced dependence on 5. In the case of "Dirac" CP conservation — <5 = and 5 = tt — the value 
of <5 is independent of the definition of the mass eigenstates (as expected), and the more cumbersome mappings are 
significantly simplified, e.g. c\ 3 -> (C23S12S13 ± C12S23) 2 / (l - c 2 13 s 2 l2 ), 5 -> 5 + it. 

It is useful to discuss a couple of well-known examples. In vacuum, the probability that a neutrino produced as a 
v a will be detected as a 1/3 is 



Pa/3 



E U <*i U l3i e 



iAml i L/2E 



13,23 /Am 2 L N 

ij=12, \ / 



13,23 



Am?,!' 



ti=12, 



(11.23) 



B [UaiU^UfrUpi] + 2 £ sin ( ) 3 [E^U^U^J . (11.24) 



It is simple, if not tedious, to check that all terms in the ij sum in Eq. (|II.24[) are invariant under all transformations 
in Table HT1 [20]. For example, a — fi, (3 = r, ij = 13 leads to 



n [u^u^u. 



t3 



l 29 



•2:i 



s 2 29 13 c 2 9 12 - c 2 e 13 s 2 e 12 



^[U^U; 3 U* rl U T3 ] = -s29 23 s29 12 s28 13 c9 13 sm5. 



1 



s29 23 c29 23 s29 12 s29 13 c9i 3 cos<5 



(11.25) 
(11.26) 



Eq. pi.26[) is, of course, proportional to the Jarlskog invariant. As discussed above, both Eq. ( pi.25[) and Eq. (III.26[) 
are invariant under 9 23 — ► 9 23 + tt or 9\ 2 — ► 9\ 2 + tt. They are also invariant under 9i 3 — > —9i 3 — tt. Finally, 9ij — > — 9ij 
accompanied by 5 — > S + tt also leaves both equations unchanged, for all ij = 12, 13, 23. 

As in the two flavor case, in order to "see" Majorana phases, phenomena involving lepton number violation are 
required. If all neutrino masses are very small (<C 100 MeV), the rate for neutrinoless double-beta decay is proportional 
to the magnitude of the ee-element of the neutrino mass matrix in the basis where the charged lepton mass matrix 
and the charged current weak interactions matrix are diagonal: 



c? 3 (c 2 2 mi + ^mje 2 * 



13,23 



s 2 3 m 3 e 2 ^ +s \ 



i ij=12 



(11.27) 
(11.28) 



It is easy to see that all terms in |m ee | 2 depend on products of s 2 9ij and c 2 8ij and are hence invariant under fty — > 9ij+ir 
or 9ij — » — 9ij. This implies that all terms are also invariant under d — > S + tt, which is indeed the case. According to 
our paraterization, all dependency on the CP-odd phases appears in the form of cos(2<5 + 203) and cos(2</>2 — 2(5 — 203). 



These also reveal that 



tt or 



tt also leave 



invariant, as advertised in Table HT1 



C. The Four-Plus Neutrino Case 



In the less familiar case of N > 1, the lepton mixing matrix isa3 + Ax3 + A unitary matrix. A general unitary 
matrix of this size contains (3 + N) 2 real parameters and at least 3 + N of those can be "field-redefined" as in the 
three generation case. Of the remaining parameters, (3 + N)(2 + N)/2 can be parameterized as mixing angles, while 
the remaining (3 + N)(2 + N)/2 are CP-odd phases. Finally, of all CP-odd phases, 2 + N will be defined as Majorana 
phases (f> 2 , . . . , <Pn+3, while the remaining (2 + N)(l + N)/2 will be denominated Dirac phases. This would be the 
end of the story if all sterile neutrino states were coupled to "sterile charged-fermions" via the charged-current weak 
interactions. This not being the case, one is free to rotate the sterile-sterile part of the lepton mixing matrix and 
further remove N 2 — N parameters, (N — 1)(N — 2)/2 CP-odd phases and N(N — l)/2 mixing angles. 

We choose to parameterize the mixing matrix as a product of all distinct R 1J (#jj), i < 3, interspersed with the 
appropriate diagonal phase matrices. Following [2l[, we define Dirac phases so that rotations between the neighboring 
planes are taken to be real, while the rest may "contain" a CP violating Dirac phase. These "complex rotations" 
in the i — j plane refer to the matrix 

R y (%,%) = P i (^/2)P-i(-^/2)R i J(%)P i (-%/2)P j (^/2). (11.29) 
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Finally, we find that we can choose all rotations in the 3 — 3 + N plane real. Using these conventions, U will be 
written as 

/3+JV \ /3+JV \ 

U 3+N = I jl R 3i (03i)R- 2i (02i, ^R 11 ^*, 6 u )j R 23 (# 23 )R 13 (#i 3 , 6)R 12 (9 12 ) f J[ P 1 ^)) , (11.30) 

where the products are defined as Yi7=i = A\A2A 3 ...A n and the expression is defined for N > 1. Several other 
choices are available (22|. Eq. (|II.30|1 singles out the "active" mixing from the "sterile" mixing and is chosen such 
that the N = 1 mixing matrix is "simple," as will be described below. 
In the case N = 1, Eq. (|II.30[) takes the form 

U< = R 34 (# 34 )R 24 (#24, ^)R 14 (#i 4 , £i)R 23 (# 23 )R 13 (#i3, £)R 12 (MP 2 (</> 2 )P 3 W> 3 )P 4 W>4), (H.31) 

^ C#l 2 c6>i 3 C#14 s9\2c9\ 3 c9\4 



where *s stand for very complicated functions of mixing angles and CP-odd factors that are not illuminating enough 
to warrant display. One will appreciate that our choice of parameterization is such that U a 4 elements, for all a = 
e, /i, r, s, and U e i elements, for all i = 1,2, 3, 4, are simple, in parallel with the parameterization of the 3x3 neutrino 
mixing matrix, Eq. pi.!5|) . The mixing matrix is parameterized by twelve real parameters: the six mixing angles 
#12, #13, 6*23, #14, #24, #34, the three Dirac phases 6, 61,62, and the three Majorana phases 4>2, </>3, 4>4- 
In the case N = 2, Eq. (|II.30[) takes the form 

U B = R 34 (# 34 )R 24 (#24, <5 24 )R 14 (#14, MR 35 (#3 5 )R 25 (#25, ^2 5 )R 15 (#15, <*1 5 ) 

x R 23 (# 23 )R 13 (#i3,<5)R 12 (0 12 )P 2 (0 2 )P 3 (03)P 4 (</'4)P 5 (05). (11.33) 

Here the mixing matrix is parameterized in terms of 18 parameters (5 2 — 5 — (2 2 — 2) = 18): the nine mixing 
angles #12, #13, #23, #14, #24, #34, #15, #25, #35, the five Dirac phases 6, 614, 624, 615, 625, and the four Majorana phases 
<^2) 03, ^4) 05- We present approximate expressions for U4 and U5 in Appendix [Cl One should keep in mind that U$ 
(and all U3+N matrices with TV > 2) is only defined up to an overall rotation of the sterile sector of the matrix. In 
the case of U5, one can always rotate around some of the mixing angles and Dirac phases via U5 — > R 45 (#45, 545)1/5. 

The next step is to define the neutrino mass eigenstates. The number of options in the N — 1 case is large, and 
we will consider only a few. The most obvious but often more cumbersome choice is to simply order the neutrino 
masses in ascending label-order: mi < 777,2 < "13 < 1TI4. A second choice — the one we will concentrate on here — is 
driven by experimental information currently available regarding active and sterile neutrino masses and mixing angles. 
Qualitatively, the data tell us that there is one state that is predominantly sterile, while the other three states are 
predominantly active. The three predominantly active states will be called u\, z/ 2 , and v 3 . These will be ordered as in 
case C3 of the previous subsection. The predominantly sterile state — by that we mean the state of maximum |?7 S i| 2 
— will be called 1/4. In more detail, |£/si| 2 , |?7s2| 2 , | C/s3 1 2 are constrained in such a way that |C/ S 4| 2 > 1/2. Translating 
this constraint in terms of mixing angles is not very illuminating sec case B3 in the three neutrino discussion — 
and will not be discussed here. Generalized versions of this can be applied to the N > 1 case, where it it is also 
necessary to specify how the "mostly sterile" states are distinguished. A simple choice is to order than in ascending 
label order, i.e., 1114 < 7715 < . . .. 

We now search for symmetries of 1/3+ jv- The procedure is entirely analogous to the one presented in the previous 
subsection, and so are the results. For all "complex rotations" R 1J (#jj, <5y), #y — > —9ij and 6ij — > 5ij + -k leave C/3+jv 
invariant, and we look for all other redundancies by looking for solutions to 

U 3+N (9' ip 6' 13 , #) = N f C/ 3 +2v(#^, S' ip #)N m) (11.34) 

where Nf , N m are 3 + N x 3 + N diagonal matrices whose elements are all possible combinations of ±1, with the 
exception of [N m ]i.i = 1. Such "sign" redefinitions of the neutrino mass states and the charged leptons allow all 



s6>i 3 c#i 4 e i5 s9 14 e lSl 
s#2 3 c#i 3 c#24 - s9 13 s9 1 4s9 2 4e- l{s+5l) c#i 4 s#24e i52 

C#23C#13C#34 - S#23C#13S#24S#34e -1152 - c6>i3s6>i4C#24S#34e~ J<51 C#i 4 c6>24S#34 



-C#23C#1 3 S#34 " S#23C#l 3 S#24C#34e- 1 ' 52 - S#l 3 S#14C#24S# 3 4e- l(,5+5l) C#i 4 C# 2 4C#34 / 

( 1 \ 

(11.32) 
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6ij — > ±(8ij + tt) (the ± sign depends on the specific value of ij) regardless of the value of all other parameters. 
Furthermore, it is easy to show that — ► —v, L identifies all fa with fa + tt. Finally, as in the three neutrino case, 
6ij — > —Oij accompanied by So- — > Sik + tt (one ik to each ij) also describes the same physics. 

In summary, we are left with the following physical parameter space for L^+at, for all values of N > 1: 

• All Majorana phases fa £ [0, tt], i = 2, . . . , 3 + N; 

• All Dirac phases 8, 8^ S [— 7r,7r[; 

• All mixing angles can be constrained to the first quadrant: cos 0y , sin 9i j > 0. 

As in the three neutrino case, we choose the Dirac phases to vary along the entire unit circle in the complex plane, 
while mixing angles only occupy at most the first quadrant. It is important to emphasize that this is a choice. One is 
also free to, for example, constrain S > 0, while allowing some of the mixing angles to also take negative values [231 ] - 

III. NONSTANDARD INTERACTIONS 

Some of the results spelled out in the previous section depend on the fact that all neutrino interactions are mediated 
by Eq. (jll.lj) . New "weaker-than-weak" neutrino interactions modify the picture painted in Sec. [ill m an interesting 
way. 

Generic non-standard neutrino interactions can be parameterized in a variety of ways. Data constrain new neutrino 
interactions to be of order or weaker than the weak interactions for neutrino processes that involve momentum transfers 
less than around one hundred GeV. Assuming that the new neutrino interactions are governed by new particles that 
weigh more than a few hundred GeV, their physics — assuming all relevant energies are low enough — is captured 
by the following effective Lagrangian, after electroweak symmetry breaking: 

£ NSI - ^ [(^V^) J» + (^a?7 Q %) A + H.C., (IH.1) 

where J~ M = A" 13 (laj^ip) + (Eaj^ep) + A^yUaj^up + A^tiaj/SysUp + ... is a generic (axial)vector current made 
up of quarks and leptons, while X — A"^(e a £p) + A"gU a U/3 + A"pU a j^Uj3 + ... is a generic (pseudo) scalar current. 
A characterizes the overall strength of the interactions (and does not concern this discussion) while £ Q/3 and i] al3 
parameterize the neutrino currents. Eq. (jlll.ip is written in the flavor basis, such that a, (3 — e, fj,, r, Si, . . . , Sjsr. We 
have also assumed that neutrinos are Majorana fermions. The rj a ^ terms violate lepton number and are further 
constrained by our current understanding of neutrino masses (see, e.g., [13] )• Constraints on for "flavor-diagonal" 
can be found in [25j ]. 

As in the previous section, we are interested in the number and range of parameters necessary to describe neutrino 
related processes. We proceed with the counting in the "standard" way. We render all charged-lepton and neutrino 
masses real and positive. We also perform as many field redefinitions are possible in order to write the charged-current 
weak interactions as in Eq. where the weak and flavor states are related via Eq. (|II.2p and U is parameterized 

as in Eq. pi.3[) . All other diagonal field redefinitions are employed in order to set U in its standard form. 

We now proceed to ask how r\ and £ change when one chooses to re-express the neutrino fields in the mass basis. If 
we assume that all r\ and £ are known in one specific basis, their values in the other basis should be properly defined. 
Furthermore, all physically distinguishable couplings should be accessible by the parameterization. 

Expressing Eq. (jlll.lj) in the mass basis, we find that £ and r\ transform as 

i -» U^U (III.2) 
r) -> U T r)U. 

In light of these redefinitions, we revisit Eq. (111.4j) . which allowed us to define equivalent classes for the mixing 
parameters 9,8,<fi. In short hand notation, U = Nf?7'N m and, of course, U^£JJ = N m ([/')tNf£Nf [/'N m . The same 
is true for r) where W is replaced by U T since both N m; f are real, diagonal matrices (Hermitian). This implies that 
any constraints lost due to £ (this will become clear momentarily) will be lost due to ry, regardless of the fact that r\ 
processes violate lepton number, so it is enough to discuss the transformation of £. 

The next step in the constraining-the-parameter-space argument is to state that 9', 5', fa and 9, 6, are equivalent 
because the N mj f matrices can be absorbed by sign redefinitions of neutrino and charged-lepton mass eigenstates. 
The new interactions impose a new constraint: 9\ 8', fa and 9, 8, <p are equivalent only if £ and Nf£Nf are equivalent, 
i.e., if Nf and £ commute. Hence, if all the elements of £ are known in one basis, the number of sign redefinitions 
that can be employed in order to constrain the lepton mixing parameter space is generically reduced. 



11 



Revisiting the discussions in Sec. [TTJ many of the parameter limitations still hold in the face of general new physics, 
while some no longer apply. Furthermore, in the N > 2 cases, one needs to revisit the sterile-sterile field redefinition 
that allowed one to reduce the parameter space by TV 2 — N parameters. We will not pursue this issue here, and will 
restrict ourselves to £ and 77 interactions that involve only active neutrinos. 

Restricting ourselves to transformations with Nf = ±1 within the 3 + N neutrino scenario, the invariances involving 
and — Oi-i t i are no longer available, while the following redundancies remain: 

1. 4>i -> 4>i + 7r,Vi > 1; 

2. 9ij — » — Oij, Sij — > Sij + tt for all "complex" i — j rotations; 

3. % — » 6 i: j + tt, 6ij + i — > -6ij+i,9 la — > -9 la ,9 bj — > -fly, for all integers b < i and i + 1 < a < j. 

Item 3 is the 3 + N neutrino equivalent of the I,R 23 (7r), — I, P 2-3 (7r) lines in Table HT1 Making use of the standard 
choices and mass-eigenstate redefinitions, we find that while we can restrict all Majorana phases to the [0, tt] interval 
and can ignore mixing angle values outside of [— ir/2, ir/2], we can't reduce the parameter space for several of the 
mixing angles any further without making choices for the elements of £ and rj.§ Ultimately, we find that we can 
choose all angles associated to complex rotations to be constrained to lie within [0,7r/2], while the remaining angles 
are constrained to lie within [— tt/2, tt/2]. We emphasize that this is a choice. As long as some angles are constrained 
to lie within [—tt/2, tt/2], others can lie in the upper-left quadrant of the unit circle. 

An example is in order. Restricting ourselves to the two flavor case, the discussion above implies that, in the 
presence of non-standard flavor changing interactions, all physically distinguishable scenarios can be parameterized if 
4> S [0,7r] and 9 E [—tt/2, tt/2], such that the sign of 9 is physically distinguishable. Further consider non-standard 
interactions where = ej^e, 2 = 0. In this case, non-standard neutrino interactions affect the coherent neutrino- 
electron scattering scattering amplitudes that govern the matter effects in neutrino oscillations. In order to compute 
neutrino oscillation probabilities, it is convenient to express these amplitudes in the neutrino mass basis: 

Aij = F z S l3 + F w U* ei U e3 + F NSI U*^U 0j , (III.3) 

where Fs are functions containing process-dependent kinematical factors and coupling constants. For our purposes, 
it is enough to know that these functions are independent of neutrino mixing parameters and couplings. The first 
two terms correspond to Z and W exchange standard model diagrams, respectively, while the third term contains the 
effects of the new interactions. In matrix form, 

F w 1 Fi 



2 





F NRI [ - (^e + 6„) S29 [(£ eAl - ^e) + {ten + £ M e) c29] 



[(£ Me - Sep) + (6m + W c29] e-*4> + s29 



(III.4) 



It is important to appreciate that only the magnitude of the different elements of ^4 maS s are relevant. This is best seen 
by analyzing the Fy/ term, which is present in the Standard Model and drives "normal" matter effects in neutrino 
oscillations. On the other hand, the relative sign of different terms within an element does matter. With this is mind, 
it is easy to see that 9 and 9 + tt describe the same physics. This is not the case of 9 and —9. This is visible in all 
flavor off-diagonal NSI terms. Also, as expected, the Majorana phase (j> plays no role in this observable. Finally, note 
that if £ were diagonal in the flavor basis, 9 and —9 would describe the same phenomenon, and the only impact of 
this NSI would be to modify the strength of the IF-mcdiatcd matter effect. 

In the case of three active neutrinos, one of the three mixing angles can be restricted to [0, tt/2] while the other two 
have an expanded parameters space. One concrete choice is as follows. For generic non-standard interactions, 

• #12,023 G [-7T/2, TT/2], 013 G [0.7T/2]; 

• 5 S [-tt,tt[; 

• 02, ^3 G [0,71-]. 



8 Diagonal but not necessarily flavor universal, implies £ = Nf£Nf , in which case the "normal" parameter space is sufficient. 
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Before concluding, we would like to emphasize that the discussion above assumes that all elements of £, including 
their signs, are known in some basis. If this is not the case, we can pick all 9^ positive and allow both signs for the 
off-diagonal elements of £. This fact was already emphasized in [12]. In the two-neutrino case discussed above, if £ 
were measured via anomalous matter effects, one would not be able to tell whether £ Me has a particular sign and 9 
is negative, or whether £ Me has the opposite sign and 6 is negative. On the other hand, the sign of all elements of £ 
should be, at least in principle, observable, as long as one can measure observables that depend on the interference 
between standard model diagrams and those involving £. 



IV. CONCLUSION 



It is well-known that while quark mixing phenomena are governed by a three-by-three unitary matrix, all physically 
distinguishable mixing scenarios are parameterized by three mixing angles constrained to lie within [0,7r/2] and one 
CP-odd phase that lies within the entire unit circle [— 7r,7r[ fl4j |. 

Assuming that only three active neutrinos exist, the situation in the lepton sector is similar but not identical. All 
physically distinguishable lepton mixing scenarios are parameterized by three mixing angles constrained to lie within 
the first quadrant, one Dirac CP-odd phase that lies within the entire unit circle 8 G [— %, 7r[, and two Majorana 
CP-odd phases that can be chosen to lie in the upper half of the unit circle, fa 6 [0,7r]. We argued that the full 
parameter range for some of the mixing angles can only be properly defined once the neutrino mass eigenstates v\, 
V2, and i>3 are unambiguously defined. In the standard definition of neutrino mass eigenstates, all three mixing angles 
are constrained to lie within Oij G [0,7r/2]. In other equally unambiguous but more cumbersome definitions (see case 
B3) smaller intervals suffice. 

In the presence of sterile neutrinos, the parameter space grows in a nontrivial way (from 6 parameters in the three 
neutrino world to (3 + A) 2 — (3 + A) — (A 2 — A) in a world with N sterile neutrinos). We find, however, that the 
physical range for the mixing parameters can be chosen in a way that mimics the standard three neutrino parameter 
space: all mixing angles can be chosen at most between [0, 7r/2] as long as all Dirac CP-odd phases are allowed to 
vary within [— tt, ir[. On the other hand, all Majorana phases fa, i — 2, . . . , 3 + A can always be constrained to lie 
within [0, tt]. 

If neutrinos participate in new "weaker-than-weak" interactions, the treatment of the neutrino mixing parameters 
needs to be reexamined. In particular we find that in the presence of new interactions which are off-diagonal in 
the flavor basis mixing angles within ]0,7r/2] and those within [— 7r/2,0[ can be interpreted as describing different 
phenomena, and the physical neutrino parameter space is doubled. It is instructive to describe a very simple two- 
flavor example to illustrate what this means. Define the mass eigenstates such that Am 2 > and assume that 
16*1 — > (7r/4)_ and that the Majorana phase vanishes.^ In this case, 

u e = (y x ± V2) , (IV.5) 

= 72 {TVl + V2) ' (IV ' 6) 

where the ambiguous sign depends on whether = 7r/4 or — 7r/4. A reasonable question to ask is whether the electron 
is the "symmetric" or the "antisymmetric" linear combination of v\ and vi. In the absence of new interactions, the 
answer is that we can't tell. Both sign choices are physically equivalent. Now add to the Lagrangian a new interaction 
of the type discussed in Sec. IIII1 where £ Me = £ e/J and all other £ values vanish in the flavor basis. In this case, the 
oscillation frequency through an electron background is modified from A = Am 2 /(2E) ( P efJi cx sin 2 (AL/2)) to 

^■matter = A 2 + + 4F 2 SI £ 2 e ± 4A(F N si£ Me ). (IV.7) 

Here, the ± sign is the same as the ± sign that appears in the electron- neutrino definition, Eq. (|IV.5|) . In this case, 
a measurement of A mat tor should reveal whether v e cx v\ + or v e cx v\ — . This statement is only true if the sign 
of i^NSi^Aie i s unambiguously determined, and the measurement of A mat ter does not do it independently In principle, 
it should be possible to measure the sign of £ Me , although it may be very challenging in practice. 



■ In the case of maximum mixing, one needs to worry about how to define the two different mass eigenstates. One way is to use the sign 
that is being discussed here, in which case the answer to the question we are asking is a matter of definition! For the discussion here, 
one can assume that \6\ is very close to 7r/4 but not identical, such that the definition of the two mass eigenstates is unambiguous. All 
results presented below are obtained under this assumption. 
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Our purpose in performing this exercise was to systematize the procedure behind the definition of the range for 
neutrino mixing parameters. Along the way, we expanded well-known (but sometimes misinterpreted or forgotten) 
results concerning the mixing angles and Dirac phase in the three active neutrino case to more than three neutrinos, 
and made explicit the origin for the [0, ir] range for the Majorana phases. We also emphasized the importance of 
properly defining the neutrino mass cigenstates before defining the physical range of the neutrino mixing parameters. 



APPENDIX A: NOTATION 



Neutrino mixing and rephasing invariances can be expressed as products of two distinct types of matrices, plus the 
identity I.* These are real orthogonal rotations in the a — b plane, R ab ((9), given by 



[R ab (0)l.. = [R ab (60] bfc = cos( 



R ab (0) = { [[R ab (0)] afc = - [R ab (0)] 6a = sin0, (A.l) 
[R ab (0)] y = %, i3¥"d> 
and single diagonal phase rotations P a (0), given by 

[pa (0) ]. i = £.. e ^. (A .2) 

Furthermore, discrete permutations of a «-> b elements of vectors are accomplished with S ab = R ab (7r/2)P a (7r). 
It is a simple matter to write an arbitrary mixing matrix with such components, as done extensively in the text 
(Eqs. (lII.5III.15III.30a 

The matrices above satisfy the following commutation relations: 

[R ab (0),R bc (0')] = sin6» (1 - cosO') A ab + sin0' (1 - cos (9) A bc + i sin sin 0' {R ac (-|) - R ac (|) } (A.3) 

[R ab (0),P a (0)] = sin0(e 10 - l) A ab , (A.4) 
[R ab (0),A ab ] = -sin(9{P b (7r) -P a (vr)} , (A.5) 

[R ab (0),A ac ] = sin0A bc + i(cos0-l){R ac (-|) -R ac (|)}, (A.6) 
[A ab , P a (0)] = \ (e« - 1) {R ab (- J) - R ab (|) } , (A.7) 
[A ab ,A ac ] = I{R b c(_|)_ R bcQ| ; (A . 8) 
while all other commutators vanish. This matrix set closes upon itself with the addition of the discrete transformation 

A ab = A ba = l R ab jpb^) _ p«( 7r )} . (A.9) 



The following identities are easy to prove and were used extensively to derive the results presented in the body of this 
paper: 



p a («/>)p a (-0) 


= I, 


(A.10) 


R ab (0)R ab (-0) 


= R ab (0)R ba (0) = I, 


(A.ll) 


R ab ((9 + 0') 


= R ab (0)R ab (0'), 


(A.12) 


P a W> + <^') 


= P a (0)P a (0'), 


(A.13) 


P»P b M 


= R ab M, 


(A.14) 


P b (7r)R ab (0) 


= R ab (-6»)P b (7r), 


(A.15) 


R ab (0)P a (</>) - P b (0)R ab (0) 


= cos0(e^- 1) A ab , 


(A.16) 


A ab P a (0) 


= P b (</>)A ab , 


(A.17) 


R ab (0)R ca (|) 


= R ca (|)R bc (0). 


(A.18) 



* One may use U(n) algebras to parameterize the n X n mixing matrix. These algebras are well-known, but their use would require 
different commutation relations for different values of n, which is not conducive to a "n-independent" description. 
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FIG. 2: P M e in vacuum as a function of #13, assuming no Dirac CP conservation. The dashed (dotted) curve corresponds 
to S = (S — 7t). If both values of 5 are part of the physical parameter space, only positive values of #13 are required to 
parameterized all distinct phenomena. On the other hand, if one also allows for negative values of 613, Dirac CP conservation 
is uniquely determined by S = 0, while 8 — n is not part of the physical parameter space. 



APPENDIX B: DIRAC CP CONSERVATION AND NEGATIVE 13 

In the three active neutrino case, all neutrino oscillation CP-violating phenomena are governed by the Dirac CP-odd 
phase S — the impact of Major ana phases is minuscule and can be trivially dismissed. CP conservation in neutrino 
oscillations would imply that e lS is real and allow a significant reduction of the lepton mixing parameter space. 
According to the standard choice of parameter ranges, real e lS corresponds to 6 = or 6 = n, such that Dirac CP- 
conservation implies that the original six-dimensional mixing parameter space splits into two disjoint five-dimensional 
ones. 

This split parameter space can be avoided by choosing a different range for 8 and #13. If one restrains 6 6 [— 7r/2, 7r/2], 
#13 S [— 7r/2, tt/2] is required to specify all physically distinguishable mixing matrices (we are assuming mass-definition 
scheme case C3) . The physical range for all other mixing parameters is unchanged. This parameterization was 
advocated in [131 ] using simalar symmetry arguments. In this case vacuum (or matter effected) neutrino oscillation 
searches can detect the sign of 0i3 US 123, or cquivalently the 6 6 [— tt/2, tt/2] versus S £ ±[7r/2,7r] value. 

Fig. fj] illustrates these two choices for the #13 parameter ranges, assuming Dirac CP conservation. It depicts P^ e 
in vacuum as a function of #13, assuming no Dirac CP conservation, for L = 300 km and E v = 400 MeV. #23, #12 and 
the neutrino mass-squared differences are fixed at values consistent with the current neutrino oscillation data [l|.The 
dashed (dotted) curve corresponds to S = (5 = tt). If one allows for negative values of #13, Dirac CP conservation 
is uniquely determined by S = — one can trivially see that the values of P^ ie associated to 6 = tt and #13 > are 
reproduced for 6 — and —#13 < 0. 



APPENDIX C: APPROXIMATE EXPRESSIONS FOR U 4 AND U 5 

We chose to parameterize the (3 + N) x (3 + N) neutrino mixing matrix according to Eq. (|II.30[) . Sec. Ill Cl contains 
an incomplete expression for C/4, Eq. (|H.32[) . Here we present approximate expressions for both U4 and U5 which 
we hope will not only serve those interested in sterile neutrino mixing but also illustrate what our parameterization 
"looks like" in the case of more than one sterile neutrino. 

We assume that all mixing angles except 623 and #12 are small. Such an approximation is in good agreement with 
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the current world neutrino data [l[. The linearized expression for U4 is 



U 3 



I 

V 



U 4 = 



-s9 12 c0 23 
s9 12 s9 23 - 




(l 

e^ 2 


V 



\ 






where 



c9i 2 

- c9 12 s9 23 8 13 e 
c8 12 c9 23 8 13 e-* s 



sVi 2 

s9i 2 s9 23 9 13 e~ 



c0\ 2 c9 23 
-c9 12 s9 23 - s9i 2 c9 23 9 13 e 



iS 
-iS 



9i 3 e tS 
s9 23 
c8 23 





( 14 e lSl \ 




9i = 


9 2 ^ 






\ ^34 j 





-cdviOue-* 1 + s8 12 c9 23 9 2i e- t52 
-s9 12 8 14 e- lS i - c8 12 c9 23 e 2i e- lS - 
-c9 23 9 3i - s9 23 9 2 ±e~ 



S9\ 2 Sf/23P34 

- c9i 2 s9 23 9 3i 



1S2 



The linearized expression for t/5 is 



U R = 



2 = 



-c9 12 8 14 e- lS ^ 
-s9 12 e u e- lS ^ 



s9 12 c9 23 9 2 4e 
c9 12 c9 34 8 



s9 23 9 2A e 




where 



723 C34 



s9 12 c8 23 9 25 e- lS ^ - s8 12 se 23 ( 
c9 12 c9 34 9 25 e- lS ^ + c9 12 s9 23 t 
su 23 9 25 e 



- lS2& - c8 23 t 
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and U 3 is given by Eq. (|C.2[) . 

It is easy to imagine how Eq. (|C.4j) generalizes to the the case of N = 3 or more sterile neutrinos. 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



'35 

'35 I 



(C.6) 



ACKNOWLEDGMENTS 



This work was inspired by a question raised by Florian Plentinger concerning the physical range for Majorana 
phases. AdG would like to thank the hospitality of the Physics Department at Columbia University, where this work 
was completed. This work is sponsored in part by the US Department of Energy Contract DE-FG02-91ER40684. 



[1] See, for example, M. C. Gonzalez-Garcia and M. Maltoni, arXiv:0704.1800 [hep-ph]; A. Strumia and F. Vissani, 
|arXiv:hep-ph/0606054| R. N. Mo hapatra et al., Rept. P rog. Phys. 70, 1757 (2007); A. de Gouvea, Mod. Phys. Lett. 
A 19, 2799 (2004J7A de Gouvea, |arXiv:hep-ph/0411274| 

[2] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 

[3] G.C. Branco, L. Lavoura and J. P. Silva, CP violation, Oxford University Press, New York (1999). 

[4] P. Minkowiski, Phys. Lett. B 67, 421 (1977); M. Gell-Mann, P. Ramond and R. Slansky in Supergravity, eds. D. Freedman 
and P. Van Niuenhuizen (North Holland, Amsterdam, 1979), p. 315; T. Yanagida in Proceedings of the Workshop on 
Unified Theory and Baryon Number in the Universe, eds. O. Sawada and A. Sugamoto (KEK, Tsukuba, Japan, 1979); 
S.L. Glashow, 1979 Cargese Lectures in Physics — Quarks and Leptons, eds. M. Levy et al. (Plenum, New York, 1980), 
p. 707. See also R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980) and 0]. 

[5] A. de Gouvea, Phys. Rev. D 72, 033005 (2005). See T. Asaka, S. Blanchet and M. Shaposhnikov, Phys. Lett. B 631, 
151 (2005); T. Asaka and M. Shaposhnikov, Phys. Lett. B 620, 17 (2005); A. de Gouvea, J. Jenkins and N. Vasudevan, 
Phys. Rev. D 75, 013003 (2007); M. Shaposhnikov, Nucl. Phys. B 763, 49 (2007); F.L. Bezrukov and M. Shaposhnikov, 



16 



Phys. Rev. D 75, 053005 (2007); A. de Gouvea, arXiv:0706.1732 [hep-ph] for recent discussions of the phenomenology of 
low-energy versions of the seesaw mechanism. 
[6] S.M. Bilenky, J. Hosek and S.T. Petcov, Phys. Lett. B 94, 495 (1980); M. Doi, T. Kotani, H. Nishiura, K. Okuda and 

E. Takasugi, Phys. Lett. B 102, 323 (1981). See also B. Kayser, Adv. Ser. Direct. High Energy Phys. 3, 334 (1989). 
[7] J. Schechter and J.W.F. Valle, Phys. Rev. D 22, 2227 (1980). 
[8] G.L. Fogli, E. Lisi and D. Montanino, Phys. Rev. D 54, 2048 (1996). 
[9] A. de Gouvea, A. Friedland and H. Murayama, Phys. Lett. B 490, 125 (2000). 
[10] M.C. Gonzalez-Garcia, M. Maltoni, C. Pena-Garay and J.W.F. Valle, Phys. Rev. D 63, 033005 (2001). 
[11] A. de Gouvea, Nucl. Instrum. Meth. A 503, 4 (2001). 

[12] N. Fornengo, M. Maltoni, R. T. Bayo and J.W.F. Valle, Phys. Rev. D 65, 013010 (2002). 

[13] D. C. Latimer and D. J. Ernst, Phys. Rev. D 71, 017301 (2005). 

[14] W. M. F. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[15] E.E. Jenkins and A.V. Manohar, Nucl. Phys. B 792, 187 (2008). 

[16] N. Haba and H. Murayama, Phys. Rev. D 63, 053010 (2001). 

[17] A. de Gouvea, B. Kayser and R.N. Mohapatra, Phys. Rev. D 67, 053004 (2003). 

[18] L.F. Li and F. Wilczek, Phys. Rev. D 25, 143 (1982). J. Schechter and J.W.F. Valle, Phys. Rev. D 23, 1666 (1981); 

J. Bernabeu and P. Pascual, Nucl. Phys. B 228, 21 (1983); P. Langacker and J. Wang, Phys. Rev. D 58, 093004 (1998). 
[19] This choice of masses was pursued, for example, in S.T. Petcov and M. Piai, Phys. Lett. B 533, 94 (2002). 
[20] J. Gluza and M. Zralek, Phys. Lett. B 517, 158 (2001). 
[21] H. Fritzsch and J. Plankl, Phys. Rev. D 35, 1732 (1987). 

[22] See, for example, D. Dooling, C. Giunti, K. Kang and CW. Kim, Phys. Rev. D 61, 073011 (2000). 
[23] M.C. Gonzalez-Garcia, M. Maltoni and C. Pena-Garay, Phys. Rev. D 64, 093001 (2001). 
[24] A. de Gouvea and J. Jenkins, Phys. Rev. D 77, 013008 (2008). 

[25] For recent discussions see, for example, J. Barranco, O.G. Miranda, C.A. Moura and J.W.F. Valle, Phys. Rev. D 73, 
113001 (2006) [arXiv:hep-ph/0512195j . A. Friedland, C. Lunardini and M. Maltoni, Phys. Rev. D 70, 111301 (2004); 
A. Friedland, C. Lunardini and C. Pena-Garay, Phys. Lett. B 594, 347 (2004); S. Davidson, C. Pena-Garay, N. Rius 
and A. Santamaria, JHEP 0303, 011 (2003); A.M. Gago, M.M. Guzzo, P.C. de Holanda, H. Nunokawa, O.L.G. Peres, 
V. Pleitez and R. Zukanovich Funchal, Phys. Rev. D 65, 073012 (2002). 

[26] D.C. Latimer and D. J. Ernst, Phys. Rev. C 71, 062501 (2005). 

[27] G.L. Fogli, E. Lisi, A. Marrone and A. Palazzo, Prog. Part. Nucl. Phys. 57, 742 (2006). 



